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VI. On the pear-shaped Figure of Equilibrium of a Rotating Mass of 

Liquid. 
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Introduction. 

This is the sequel to a previous paper on “ Ellipsoidal Harmonic Analysis ” (‘ Phil. 
Trans..,’ A, vol. 197, pp. 461-557). I here make use of the methods of that paper, 
and for brevity shall refer to it as “ Harmonics.” 

The sections 1 to 4 are preparatory, and might have been included in “ Harmonics,” 
but seem more appropriate here. Section 5 is an independent investigation of so 
much of M. Poincare’s celebrated memoir on rotating liquid* as relates to the 
immediate object in view. 

It is not necessary to say more here, since 1 give a short summary in the last 
section. 

§ 1. The Harmonies of the Third Degree. 

It was remarked in § 7 of “ Harmonics ” that this group of harmonics may be 
determined in rigorous algebraic form ; I now proceed to do this. 

The requisite formulae are contained in § 6 of “ Harmonics,” and I take the several 
harmonics in succession. Throughout this section we have, of course, i — 3. 

s = 0; type OEC, and $ s (v) = q 0 P s (v) + with q 0 = 1. 

The equation for /3a is 

_ jfifl} {3,2} 15/3 2 

1 “ 4 + /3a 1 4- p3a ’ 

If we write 

{B x f = 1 +15y8 3 , 

the proper solution of the quadratic equation is 

/3a = 2 (B x - 1). 

* 1 Acta Mathematical vol. 7, 1885. 
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Then since 
we find 


TXfc = 2&Trj’ with *» = L 

•§P — fj/, P.f 2 = 15v (P — 1), 
2/9 — 1 + — 5 ( 1 — 13 )) 


= 1; type OOO, end P» = («iW (-) + l»*W M), with ft' = !• 

The equation for /3<r is 

*r + U>. 3.4 = 


If we write 


' 1 8'+ /3ff 

_u> pa 

or . Ar+6/5 = r+ i^. 

If we write 

(&)*=i-f/8(i-/8), 

the proper solution of the quadratic equation is 

0<r = 4(3*-1-10. 


?£s. — . . if:__J:.. ai1 q ?s_o £s with o' — 1 

?1 ~ 8 + f3a ~ 4-(P, + 1 - f/8)’ d ft' ” d ft’ rtl _ • 


(5P - 

- 1) (P — l)' 1 , 

P 3 3 (?) = 15 (P - 1) (P - 1 

.)*, we find 

(*»- 

- 1 + 2/8) (P 

1 + /9\*/o 3-/3-2RP 
i — ^9/ V Td-ft’j 

. . .(2), 


.s — 1, type OOS and iPg 1 (v) = tftPo 1 (v) + (3q :] P 3 S ( v ), with q x = 1. 
The equation for /3cr is 

o_ ,o o d _»[M 


/3cr — J/3 .3.4 = 


8 + /So- 


A comparison with the last case shows that we have only to change the sign of /3 
accordingly if 

(l\f =l+f0(l + j8), 


we have 


/3cr=4 (Jig - 1 + fj 8 ), 


ft" 4(5“+ l+f/8)’ 


with q l — 1. 
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On substitution we find 


»,» = J CB, - 1 + 3/3) (v> - 1 )*(k - bf( 


» = 2 ; type OEC, f 3 3 (y) = fiq 0 I\ (y) + q 2 P^ (v), with q % = 1. 
The equation for /3cr is 

„ -i/8»{3,l}{3,2} _ - 15/3 2 


/3cr = 


4 —/3<r 


1 — f/3 <t ' 


We have already defined (Bff as 1 + 15/3 3 , and find the proper solution of the 
quadratic equation to be 

/3a = —2(^-1). 


%„ _ - {3, 1} {3, 2} _ - 4 (% - 1) 


4 -/3<r 


-, with q% = 1. 


With the known values of P> and of PJ, we find 


®» , « = v( 1 + S ^-- B i)v(v*- 6 4 ^) .(4). 


A comparison with (l) for s — 0 shows that the last factors in each only differ in 
the sign of B v 


s- 2 ; type OES, P 3 2 (v) ~ V • P a 2 H 

Since id, 2 (v) = 15y (y 3 — 1), we have at once 


P 3 3 (y) = 15y(y 3 - if 


,3 1 ± ft 

1 -/ 3 , 


/ a I+ft 

3 ; type OOC, P 3 3 (y) = a/ [fiq.'Pfiv) + q. A 'P 3 s (y)], with q.' 

V* — 1 


The equation for /3cr is 


We have already defined 


W{ 3, 2} {3, 3} 

8 — /3<r — 6/3 8 — /3cr — 6/3 


a ) ! = 1 - i /8 (1 - / 8 ). 


and find for the proper solution 


/3cr =4(1 —-1^8 — i? a ). 
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Also 


and 


2ft ~ {3, 2} {3, 3} 
ft 8 - pa - 6/3 




£l 

ft' 


5'-, with g./ 
3 ft 


1. 


Whence on substitution 


p 3 8 M 


20 

0 


(1-5.) 


1. + /8V / , 

i-/s) r 


,2 


3 — /3 + 22? a \ 

5 (1 - >) / • 


( 6 ). 


• s ' — 3 ; OOS, ^ 3 3 (v) = 0 q x Pg 1 (v) + g 8 P 3 3 (v), with g s = 1. 
The equation for yScr is 


/3<X 


i/3»{3, 2} (3, 3} 


8 - /3<r + 6/3 ' 

We may derive the result from the last case by introducing 

(P3) 2 = 1 + f /3 (1 + 0 ), 

and changing the sign of /3, so that 

y8cr = 4 (1 + §0 — Pg), 

2 A _ 16 


ft 


Whence on substitution 


/3 2 


(P s — 1 - i/3), with q 3 = 1. 


1 (,# , o + <8 + 2B 3 


®. s w= (* + f> - b,) p - a* y 5 (1 _ fi) 


(?)• 


and 


The forms of the corresponding functions of p are the same, except that (1 — p 3 )* 

1 + j3 \i 

l~ZTp — /P) replace the corresponding factors. 


I have not determined the cosine- and sine-functions, because they may be written 
down at once from the results already obtained. The three roots of the fundamental 

cubic are zd, pd, and - . Hence we have only to replace zd by this last 

function in the seven fbrmulse (1)—(7) in order to obtain functions proportional to 
the seven cosine- and sine-functions. If the definition of the latter functions is to 
agree with that given in “ Harmonics,” the factors must be determined appropriately, 
but the question as to the value of the factor will not arise here. 
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§ 2. Change of Notation. 


It will be convenient, with a view to future work, to change the notation, and I 
desire to adopt a notation which shall not only agree in the main with that used 
in “ Harmonics,” but shall also facilitate reference to a previous paper on Jacobi’s 
ellipsoid (‘Roy. Soc. Proc.,’vol. 41, pp. 319-336). 

I write 


_l_-0 

1 + /9 ’ 


Ji — 


It may he noted that what I here write k was denoted by k in “ Harmonics, ’ and 
vice verst 1 

I have in general written the current co-ordinates v, g, <f>, and the ellipsoid of 
reference r 0 , so that the squares of the semi-axes are 



1 + P \ 
i -fir 


‘( v 0 2 ~ l )> 


I now propose to write as the squares of three semi-axes of the ellipsoid of 

y»p-pPT*P'n PP 

c 2 cos® y, c 2 (l -K 2 sin 2 y), c®. 


Comparing these two we see that 

k = ck sin y, and v n = —r—. 

' u k Sin y 


For the current co-ordinates I retain and write 


v 


1 

k sin ^ ’ 


ji — sin 6. 


The three roots of the fundamental cubic are therefore 


A 


sin 2 6 , —\ (1 - k* cos 3 <f>). 


/r sin 3 ^ ? f ~ .‘ 1-/3 

The rectangular co-ordinates x, y, z are therefore now expressible as follows 


x == . cos xb (1 — k 3 sin 2 0f cos <A, 

sm x _ 1 

c sm 7 f /i • » 

v = 7—- . (1 — snr i by 1 cos # sm <£>, ) . 

° Sill 

2 = C’dTT _ g ; n q /•) _ K 't cog 3 ii} 

Sill -yfr ' T/ 

2 R 


(8). 
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30(5 
These give 

af , _ f __ _i_, 3 _ r L Kil|J 7 

cos 3 ip 1 — k 1 siu 3 \j/ Z siir 

At the surface \fj = y, and we have 

x* » 2 , „ 

a h T"~i" 4“ 2 — c". 

cos 3 7 1 — sin- 7 


In the formulae for the third harmonics, in every case hut one, and in two out of 
the five harmonics of the second degree, there occurs a factor of the form (iv 3 —■ con¬ 
stant) ; in each such case I write that constant in the form ^ 3 /k 3 , and q' z = 1 — <f. 
Thus q will have a different value for each harmonic. 

It has been already remarked that for most purposes it is immaterial by what 
constants the several functions are multiplied. Although it would be easy to 
determine the constant in each case so as to make the function agree with its value as 
defined in “ Harmonics,” yet I shall not take that course, and shall omit factors as 
being in most cases redundant. 

For the sake of completeness I will give the first and second harmonics in the new 
notation, as well as the third. 

Since the harmonics of the first degree are expressed by 

f i (v) = v , Pji W = 

it is clear that in the new notation 

= P ‘*W 

f i (p) = sin 0 , Pi 1 (fi) 

C 1 (<^) = (1 -k'W*)*,' Wit) 


r|)‘. 

= c ot x h P/ (v) 

= (1-kW 0)*, ^(p) 

= COS < f ), S] 1 ( <f>) 


_ (1 — n 2 sin 3 | 

sin yfr ? I 

— cos 6 , K 9 )’ 

= sin (j f) ) 


It appears from § 12 of “ Harmonics ” that 


is (v) = V 2 + 


7 

3 

a 


W 0) =+ 


/ 

7 

/ 3 
Ot 


where 


7 It — 2 y' 

a ~ !*, (1 - £) ’ «' 


5 + 2 
:>(! -£) ’ 


and 


ZF - 1 + 3 y 8 *. 


In accordance with the notation suggested above, let 
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# « 1 —• 

Then substituting - — for /3, we find 

-I. ~f~ /C~ 


and for both cases 

Hence 




9_ 

9 9 " 

k 1 = q z -- 


Pa (*0 and $ a * (v) = 


1 — 2<? 3 ' 

1 — (f sill 3 yjr 
sin 2 yjr 


^ 2 (/x) and $/(/*) = ^ 1 ~ 3 sin 3 6 ), 

and C/ (</>) = 1 - *, 3 cos 2 </> 

where «* — r / 3 |-- , and <f — - 3 - [I + k 2 (1 — kV 3 )*], with 

first and the lower sign for the second. 


It appears from (19) and ( 20 ), § 7, of “ Harmonics” that 

P* 1 (*0 = ^77, 

" v ' Sill 3 yjr 

Pg 1 (jx) = sin 6 (1 — /c 3 sin 2 Of, 

Cs 1 (<f>) = COS (j) (1 — k ' 3 COS 3 (j)f 

and from ( 21 ) and ( 22 ) that 


f 3 ‘ w = 


(1 — /t a sill 2 yfrf 


sill 3 f 
Pj 1 (/j,) = sin 6 cos 0 , 

S. 3 1 (<f>) = sin (j) (l — /d 3 cos 3 (j)f 


from (25) and (26) 


P.*w 


00 s \jr (1 — /c 3 sill 3 yjr)* 


sin 3 yjr 

Pf (jx) = cos 0(1 - k 3 sin 3 0) a , 
<Sy (r/>) = sin (f) cos <f> 


Turning to the harmonics of the third degree, we found that in 


the type is OEO, 


P i (") and p 3 2 (v) — v( v- 


4 + i>\ 
5(1 


%) 


2 it 2 


• ( 10 ), 


upper sign for the 


• ( 11 ), 


• ( 12 ). 


. (13). 


the* two cases Avhere 
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If we put 


t _ ±±A p 

* s ~ 5 (1-/3)’ AVe " H 

( / = f[l + K 3 T(l-fK 2 +K 4 ) 5 ], 


« , 4 - 5 r/ 

K — (7 , - -— 

1 3 — 4<f 


Therefore, with the above alternative form for cf, 

f 3 (*) and W (v) ~ i "' r ’ ,S ' U ^ * 


sin 3 ^ 


(/x) and ( i a) = — sin 0(1 — — sin 5 0 j > 

C ;! (</>) and C/ (<f>) = (l - ” cos 3 c/>) (1 — /c 3 cos 3 </>)= 


Again in the two cases where the type is OOC we found 


0 8 — /3 4 - 2A 3 


V O') and P 8 »(v)= k 


5 (1 — /3) /V 1-/3; 


Putting 


(/ _ 3 - /3 + 2/1, 

« 2 ” ' 5 (1-/3) 5 

2 3 = -H1+2k 3 T(I-~k 3 +4k^) 


9 ,2-5^ 

K= ^' 


Therefore, with the above alternative form for f/, 


P; (v) and P t > W 


cos ig (I — sin 2 yjr) 
sin 3 A 


P ;! l (,x) and Pg (g) = - (1 - k 2 sin 3 6f [l - - sin 3 6 


C3 1 (<£) and C3 3 (</>) = cos <f> ( 1 - ^ cos 3 c/> 


In the two cases where the type is OOS we found 


f ; (a) and f ; (a) = (b - “rljf ) (» S - 1 ) ! 
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Putting 


we find 
and 


3 + 0 + 2 B, 


/c 3 5 (1 - 0) 

^ = |(2 + ^qF( 4 - K V 2 )*) 


k* = q 


,4-5^ 


1-2 cf 

fherefore, with the above alternative form for q 2 , 


« -d «,> (,) = , 1 


W (/*) and $ 3 3 0 *) = - cos 0 ( 1 — ^ sin 0 6), 


S ;; 1 (</->) and S 3 3 = sin <f> (1 — -- cos 0 <J) 




(16). 


The seventh of these harmonics, which is of type OES, stands by itself. We 
had 

P/M = -(v 3 -l )*A - 1 + W 


1 -0,1 


This gives in the new notation 


p 2 /„\ _ cos-f (1 - , 

3 ' ' sin 3 -^ 

P 3 0 (jx) = sin 6 cos 6 (1 — k 1 sin 0 0)- , 

S 3' 3 (<f>) = sin (b cos <fy (I — k 3 cos 2 (/>) J 

The formulae (9) to (17) give the fifteen sets of three functions constituting the 
fifteen harmonic functions of the first three degrees. It would be easy, although 
somewhat tedious, to find the coefficient by which each function is to be multiplied so 
that its definition may agree with that of the previous paper. 



§ 3. Expressions for the Solid Harmonics in Rectangular Co-ordinates. 

The three roots of the original cubic equation were v 2 , /x 2 , — - , and in the 

new notation the three roots of 


r 


I ,7 I v> O * .7 1 - 7 .1 1 K ' J COS- (j) 

— — TJ-, + V , + — = cVsur y are - , snr V, -— 

w 3 — 1 /k 1 or — 1 or rsm-f k- 
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Hence it follows that we have the identity 

« 2 1 / 2 ~ 0 1) • 0 .7 • 3 (tirv — co") (sill." 6 — oj-) (i-.-AWi — w -) 

- -,-p, +. 0 , + .1 - sop y = C-K- S1I.P y -WL ... 

w- — 1 Ik~ or — 1 ft)~ (1 /«■' — w") (1 — ft>")ft)" 

Putting co 1 = %,, 

2 3 . 0 . „ C 3 sill 3 7 , 0 • o ,W, i * 3 . 0 /, Wt *' 2 0 j\ 

4. A __ - -f c"snPy= r - 0 — 3 ... (1 — (f sup i//)(l — -,sm 3 0) (i — 75 cos~ <jf>). 

c / 3 k : 1 — q~ r/ ' (k~ — (f ) sni" v J - ,/v 5 " /v r/ 

This expression, together with those for *, y, z in (8), enables us to write down the 
results at once. As before, I drop the several factors as being redundant for most 
purposes. 

From ( 9 ) 

IMOfiMOR) = 3. P, 1 M P,‘ Mffi, 1 «■) = *, IVMS, 1 MS, 1 «,) = >/ 

... ( 18 ). 

From (10) 

• S. W IPs M «. W »»d ip. 3 M IPs 3 w «, 3 w = jw + »• - *-%• 


4 - c 3 gY a s in 3 y . ( 19 ), 


(/ = l[l + k- ^ (1 - kV 3 )*] , and k 3 = 7 


1 - 


so that 


; 7 = 1 - 

K‘ — q- 


From ( 11 ), ( 12 ), and (13) 

P, 1 (p) p, 1 00 C4 (*) = a*, w (p)-W 00 s, 1 (+) = yz , Py (0 Py 00 (+) = xy 

• . . ( 20 ), 

From (14) 

$ 3 (0 fa OO 0 3 (</>) and f 3 2 (*) f 3 a 0*) C 3 ~ (</■») = z (<i % 2 + r ~ tf 2 * 2 

+ (r<fq~ sin' y) . ( 21 ), 


where 7 3 = f [1 + k 3 4- (1 *— i k 3 + #c 4 ) 4 3 » and k 3 — 7 ".! _ ^ 


so that 
From (15) 


“Ts = 3 — 4 q\ 

K* ~~ Cf 


Pa 1 (p) P a 1 00 C, 1 (*) and P a » („) P 3 * (p.) <£/’ (+) = * (JV + ^ ?/ - r 


-)- sin 3 y) . (22), 
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where 

so that 
From (16) 


<f = |(1 + 2 * 3 + ( I - + + 4k 1 )') , and k % = <f 


a 2 - 5 r 


1 - 4g» ’ 


T= 1 - 4g 3 

K" — ' t 


(f(p 


»' W Ws 1 0*) s, 1 (<» and ft 5 (-) *>* M S,‘W = y (l V + / - «'V 

4- c 3 gY 3 sin 2 y) . (23), 

, 4 — 5§ 3 


where i/ 2 = 4 (2 + k 2 4- (4 -- kV 3 ) 5 ) and k~ = 5 a -j—^ 

so that ' "l£--- 3 = g (1 — 2(f ). 

Lastly, from (17), 

P 5 i WP.*W8 s , W = ^ . • • • 


(24). 


It is easy to verify that each of these expressions satisfies Laplace’s equation. 


§ 4. The Expression for the Q-fmctions in Elliptic Integrals. 

In this paper I drop the factors (£ and E which were found to be necessary when 
the Q-functions were expressed in series. 

We make the following definition :— 


Wiiy 0 ) ©f+o) — [f/(p 0 )] 3 f 

j V, 


civ 


[$/(+] 


and a similar formula holds for P/ Q •. 

It is clear that may be multiplied by any constant factor without changing 
the result; hence we may use the forms which have been found in §§ 2, 3. 

The notation must now be changed. 

I 1 

We have v — ——, and n n — —:— . Therefore, when i h is adopted as variable, 

K Sill Y K 81117 T 1 ’ 

the limits are y to 0, and the sign of the whole is changed. 

But 


and 

Therefore 


O' 3 — 0 ; \ v 


J cos 4 7 , 

civ = — . , cl\p, 

K Sin- yfr ' 

1 + fi\i _ cos 4 (1 — k" sin 3 4)- 


3 - A 


ic~ sin 2 4 
d\Jr 


r dv _ [V cl- 

hv~- D* +-r~P' “ * ] 0 (i - k 2 


; sin 3 a|t)' 
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In accordance with the usage in elliptic integrals, I write 

A 3 = 1 — k 3 sin 3 xjj 

under the integral sign, or 1 — k 3 sin 3 y outside the integral. 

I shall also for brevity write 

A t 3 = 1 — (f sin 3 tf/ 

under the integral, or 1 — q 2 sin 2 y outside the integral. 

We have then 

»<• w ®? w=«»• w? • 


I apply this formula successively to the several functions, as given in (9) to (17), 
and introduce the abridged notation just defined, but I do not reiterate the special 
meanings to be attached to the symbol q in each case. 

Since Po( p ) ^ 1) we have (dropping the now unnecessary suffix 0), 


PoWSoW = . 

Pi , WQ, 1 (-) = «or r f’^ > |±#, 
Wr) ®< M <»d V S W «,* (-) = }’ 


_ k cos 3 7 A 3 p sin 4 'fr 


P*(v) Q*tr) 
z \ / \ / sin* 7 j 0 


J A COS" y!rA a 7 


P^WQs 1 


P s s (-) Q,»W 


yJrA 

kA l 1 p s.in° yfr 
sin 6 7 J 0 A 1 ‘A ( A i f ’ 

~ a:cos 3 7A 1 4 p sin 6 \jr 


*y 
* n 


sm f) 7 J 0 cos" 

K « ® s ‘M «>d ».•<-) « S 8 M 1' gj #. 


p 5 *we. 


/c cos 3 yA 3 [v sin 6 


sin 6 7 


[y sin 6 

JoCos^A"^ 
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All these integrals are expressible in terms of the elliptic integrals 

m:?. ■*=]>* n =- 

It will, however, be found that in fact the coefficient of II vanishes in every case. 
The cases of i = 0 and i = 1 are very simple, and we have 

Mo = 

(>->), 

P^Qi 1 = x cot 3 y (\ A tan y — ~E), 

\ fC fC j 


30 1©/ = a ' A! (- v F + -±- E - Hi " 7CnS7X i 

Wl ^ sin 2 y \ « 2 ^*V 2 *'*A / 

It is possible by direct differentiation to verify the following results, although the 
verification will be found pretty tedious. 


f »m 4 t rl , _ (2 - - Y(1 - 2f) 2- 1 

J A X 4 A V — 2gW - g 2 ) "*■ 2 2 Y 3 ' 

f sin 4 1 1 — 2/c 3 T , , 1 

I cos 2 ^A C ^ ~~ ~ Y 1 + VW E + Y 3 A tan ^ 


1 p A sin yfr cos t/t 

27w-2 3 ) + ’ 


«V 3 A 


fsin 4 ^ 7 2 , 1 + /c ' 3 ^ sin \lr cos \lr 

j-*? #=-?* + - 

f sin 4 -^ , _ 1 l + /c 3 tan\|r 

J cos^AS ^ “ /eV 3 ^ M + *' 4 A 


•2^'2A ’ 


*=«* * - y?- * + ;.r [*-(»+«*> »■« 

These are all the integrals needed for the harmonics of the second degree. In the 
case of the first we have 

Thus the coefficient of n vanishes and the results are 

F- l ~^E + 

P s ’ W Qi « = [- «. F + Wf ^ ^ A bmy J , 

™ l /.a ^ i / \ _ * A3 f 2 v , 1 + «' 3 V sin y cos yl 
$2 (0 ©3 (0 - sin 4 y L Y * + «V 2 Jj ~ ^VAA~J ’ 

P 3 It\ n 3 hr\ K cos27A2 1 J? ^ "h ^ T/ 1 1 * an 7(2 — (1 + Y) sin 2 y) 

^2 Wy» 10 - sin 4 y I A ' 21 ~ A + AA * 


In the first of these <f = -|[l + Y ^ (1 — YY 2 ) 4 ] and Y = <? 3 


2-3g » 
1 - 2 ? 2 


VOL. OXCVTir.—A. 
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The following integrals may also be verified by differentiation : 


si B!± ,1, = ti 4 - m - * 3 (3 - 4g 3 ) 2g 3 g ' 3 + A 2 (3 - 4 g 3 ) __ zc 3 (3 - 2g» ) - 2g 3 g' 3 

Aj 4 A • 2gV 3 (zc 3 -g 3 ) " 1 ” 2zc 3 gY 3 2zc 3 gY 3 (zc 3 - g 3 ) J 

A sin cos \[r 


sin 8 ijr 
cos 3 'i|rA 1 l A 


dxp 


^ 2g 3 g' 3 (zc 3 - g 3 ) A/ ‘ 

zc 3 (l - 4g 3 ) - g 3 (2 - 5g 3 ) 4g 3 - 1 1 + 2g 4 - zc 3 (2g 3 + 1) 

2gV‘ (zc 3 - (f) 1 + 2gY 4 "* ~ 2 2 ) 

A tan i/z- A sin cos i|r 
+ "zdy 4 “ 2g /4 (zc 3 — g 3 ) Aj 3 ' 


(26), 




(27), 


f «hdjf , __ zc 3 (1 — 2 g 3 ) — g 3 (4 — 5g 3 ) 2g V 3 - zr (l - 2g 3 ) _ 2g 3 g' 3 H-jAe'® 

] A X 4 A 3 2gY 3 (zc 2 -g 3 ) 3 2zc 3 gY 3 (« 3 - g 3 ) 2 zc 3 zc' 3 g 3 g' 3 (zc 3 - g 3 ) 3 ^ 


sin ^ cos i]s A sin i/c cos ap 

^zcT^Tg 3 ) 2 A “ 2g' 2 


sin 6 \js 


d\jj = 


2 — ir p 2(1 — zc 3 zc' 3 ) v ( sin i/c cos i// ( A tan 


J cos 3 A 3 

Now in (26) we have to pnt 

in (27) 

and in (28) 


zc 4 zc' 4 ^ zc 2 ze' 4 A 




zc 2 = 


= 2~ 


.s 4 "^ 2 . 
3 - 4g 3 ’ 

2 - 5g 3 _ 

mg 3 ’ 


k==y- 5?! 


1 - 2g 2 


(28), 

(29). 


Introducing these values, and taking the integrals between the limits y and 0, we 
find: 


33 q ;U1C [ 33 any 2 -«A, 4 J7g' 2 1 p _ 2g 4 + og'-_1 ^ (4gy-1) A_sin 7 C 0 S 7 | 

IfWs ancl IPs <9-3 — a ,-,,6 .,1 oTwU V 9.^rAn'i ' 2g 4 g' 4 A 3 I 


P^Q^and P 3 3 Q3 3 = 


sin 8 7 1_ 2ze 3 g 4 g' 3 * 2zc 3 gY 4 

zc cos 2 7Adf 4g 3 - 1 r 1 - 5 g 3 - 2g 4 


Ad J 


(30). 


sin 0 7 [ 2g 4 g' 4 




V4 


J? 


1 — 7g 3 — (1 — 5<2 3 — 2g 4 ) 


Sn'i 


^ sin 3 7 ~j A ta n 7 | 


ft 1 ®,' and ¥, 3 ®, s = ^ J + E 


y 3 

6zcV 3 gY 

1 — 5g 2 + 6g 4 - g 3 ( 2 — llg 3 g' 3 ) sin 3 <y \ sin 7 cos 7 ] 


6 zc /3 g 4 g ,/4 


Z \ 0X11 7 7 ! /oo\ 

V^Ui ■ (32) - 
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p 3 2 Q 3 


k cos 3 7 A 3 f 1 + K' 


sin 6 7 


F 


i(l 


E 


+ 


1 4 - $ — (1 + kF) sin 3 7\ tan 7] 




} • (33). 


In (30) 

q* = f[l + k» =F (1 - {k 2 + K 4 )*], 

2 3 - 4 2 3 

In (31) 

(f — i[l + 2 k 2 “F (1 — k 3 + 4/c 4 )*], 

9 9 2 - 5 s» 

K =5 3 i _ 4?3 

In (32) 

f = i [2 + (4 - kV*)*], 

2 0 4 

K ' = 


§ 5. Bifurcation of Jacobi’s Ellipsoid. 

If a mass of liquid be rotating like a rigid body about an axis, x, with uniform 
angular velocity co, the determination of the figure of equilibrium may be treated as a 
statical problem, if the mass be subjected to a potential -|co V + s 3 ). 

The energy lost in the concentration of a body from a condition of infinite disper¬ 
sion is equal to the potential of the body in its final configuration at the position of 
each molecule, multiplied by the mass of the molecule and summed throughout 
the body. In the proposed system, as rendered a statical one, it is necessary to add 
tt&t (if -j- z 2 ) to the gravitation potential before making the summation. If A denotes 
the moment of inertia of the body about x, this latter portion of the sum is |-A co 2 , and 
is therefore the kinetic energy of the system. 

If dm l5 denote any pair of molecules and D n the distance between them, and 
E the energy lost, we have 

E = + iAm > 

J Ui 2 

If the system had been considered as a dynamical one, the expression for the 
energy of the system, say U, would have resembled that for E, but the former of 
these terms would have presented itself with a negative sign. 

It is clear that the variation of jrAco 2 , when the moment of momentum is kept 
constant, is equal and opposite to the variation of the same function when the 
angular velocity is kept constant. 

The condition for a figure of equilibrium is that U shall be stationary for constant 
moment of momentum, or E stationary for constant co, in both cases subject to the 
condition of constancy of volume. The variations in question lead to identical results, 
and I shall proceed from the variation of E. 

2 s 2 
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Jo (u + a?f (u + Wf (u + '-O’- ’ 

the internal potential of an ellipsoid of mass M and semi-axes a, b, c is 


I ** + £** + **?- 

1| T «(& T } clb ^ c do r 


Hence 


dm, dm* 


Qfffb i ( 

^ j i) 


»<[*+?£+ • • •]<*»• 


Now if A, B, C denote the principal moments of inertia of the ellipsoid about 


x, y, z, 


j* x\lm = %(C + B -A) = lMa\ 


and similar formulae hold for the two other axes. 
Therefore 


dm.dm.xy 07,^0 [% , if c ^ r . 1 ^ . c ^ r 

*] -ir iU " r + * r&+ 6 * + c *. 


But since 'P is a homogeneous function of degree — 1 in a, b, c, the sum of the 
three differential terms is equal to — TC Hence this expression is equal to y () MA'. 
Since 

iAco 2 = -^M(b z + c z )co z , 


we have 


JS=&AP + 


If E be varied, whilst abc and w are constant, it is stationary if 

<W s . ftW . 2b A /dW .2c \ 

— Za + + — *>~j 8b + + — co ) Sc = 0, 

8a 8b Sc _ 

“T r ' "T — 0< 

abc 

Eliminating 8a, 86, 8c we have the well-known conditions for Jacobi’s ellipsoid 


2« 2 6 3 _ <AP _ j dV 

3M a cla db’ 

2 ®*c* _ d9 <AP 

3AI ~ Ci da ~~ C d? 

1 / fhP , fMA 1 / <AP AP 
nr a —— 6 — = — a — — c t 


• (34) 


A 3 v da da 


\ da ° dc I ' 
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If we add together the first two of these, and avail ourselves of the property that 
9 is homogeneous of degree — 1, we easily prove that the stationary value of E is 


E = if 3 j~T -f- a 


Since the potential of the ellipsoid must satisfy Poisson’s equation 

09 09 09 _ 2 

ada bdb cdc abc ’ 

Also 

, 09 . 09 

a E, + h E +C Tc~- ,qf - 

By means of these and two out of the three equations (34), we may eliminate the 
differentials of \P, and writing p for the density find 


+ *- + j) - 6 

(&2 + C 2 ) + p- + ~) ~ f 


■ • • ■ (35). 


I do not happen to have seen this form for the angular velocity of Jacobi’s ellipsoid 
in any book. 

It is easy also to show that the stationary value of E may be written 


E = -i s M* 


1 . 1 


( ( - ! + cS )( i +^+^)- 6 


¥ + c 2 


We may now express the potential, say V, of the system entirely in terms of 9 

and a , for 
da 


V — \M *+^ + Sr« 


a? d9 y l ( 09 2(oW\ z* / d9 2®»e* 


a da ' ¥\ da 3 M 1 ¥\ da m 


4 * "5 ( a 


+ W (y' 2 + z ~) 


■iuW + a^d + t + i, 


We thus verify that V is constant over the surface of the ellipsoid. 

Let g denote the value of gravity at the surface. Then if dn be an element of the 


outward normal, g = 


Since 


dx px dy py dz pz 

dn a 2 ? dn b 2 ? dn c 2 5 
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where 


- + • / -+ - 
,A ' U ~ A ’ 


9 


, %r dw fa* f <P\ 

° Ma ’ + w + ? 


da 


, . M dW 

2 .« ",' 

p dci 


Now change the notation and write 


Then 


Now 


and 


«? = V (v - j 4 |) > V = P(-o* - 1). 0 s = *V, 

M = /c 3 (v 2 — V 0 2 ) . 

2 r _ <&,_ 

_ * - rr^D’O 3 - !)* 

^ _ 2 / 8 1 + /3\ f cfc 

® d*' ” * v° i - J ,„(" 2 - !: ! J) 5 <v ! -1)*' 

Po W = 1 » Pi 1 W = (** ~ ’ 

f oe 

F [fTppTIT-p) 5 yiTriji ; 


so that 


and 


^ = f $0 ( v o) ©0 K) y 

»f=- t 2p , , WQ,'W, 
g = l f Pi 1 M Q, 1 « J 


( 36 ). 


We may note in passing that the condition for a Jacobian ellipsoid (the last 
equation of (34)) is reducible to the form 


«A 2 

sin 4 <y 


dip = K cot 3 y r 
J A° J, 


tan 3 -^ 


dxji. 


On examining the series of functions given in (25), we see that it may be written 

W ( v o) ©2 1 ( v o) = Pi 1 ( v o) Ql 1 ( v )o • 

This agrees with M. Poincare’s equation (l) on p. 341 of his memoir. 

We will now suppose that the body, instead of being an ellipsoid, is an ellipsoidal 
harmonic deformation of an ellipsoid, which is itself a figure of equilibrium for 
rotation eo. 

The addition to E will consist of three parts; first that due to the mutual 
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energy of the layer of deformation ; secondly that due to the ellipsoid and the 
layer ; thirdly that due to the change in the moment of inertia. 

If a subscript l denotes integration throughout the space occupied by the layer, 
U the potential of the ellipsoid, and dv an element of volume, 

SE = + [ Updv + \(f + 1?) pdv. 

If £ denotes the thickness of the layer standing on the element da, the first 
of these terms is ^ p 3 • 

The value of U + (y 2 + z 3 ) throughout the layer is equal to F 0 — gt,', where 

V 0 is the constant value of U + 4 « 2 (y 2 + z 2 ) over the surface of the ellipsoid, and 

is the distance measured along the normal to the element d'Qda of volume. 


Hence 


j ' Updv+ Jw 2 J ' (y 2 + z*)pdv = Jj p(V 0 — gQdl'da. 


Since V 0 is constant and the total mass of the layer is zero, this is equal to 


T" j 'g{? da. 


■ fiP 
It follows that 




The axes of the ellipsoid have been chosen so as to make our original E 
stationary, and the further condition to be satisfied is that 8 E shall he stationary. 

Let us suppose that 

C = peW(p)^i s (</>), 

which expression shall be deemed to include any one of the other types of harmonic. 

Then it is shown in (5L) of “Harmonics” that the potential of this layer at the 
surface of the ellipsoid is 

■? 1/ 

e (v 0 ) €14 Oq) $/ (/x) <BV (<£). 


Since the mass of an element is pep (p) (&f (cj>) do -, we have 


1 

2P 


^I ^3^^" ~\d ^ 9 

D~ 


3 . 

2 k 


Me.* 


W ( p o) (uj) f [W (p) ($)~f pda. 


With the value of g found in (36) 




} gedo- = |ye’P, 1 WQi 1 (*o)f [%’ W <M (■»]* pdcr. 
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Hence 

In order that the new figure may be one of equilibrium, this expression must oe 
stationary for variations of e. It follows that we must either have e = 0, which 
leads hack to Jacobi’s ellipsoid, or else 

i _ (bo) On) 0 

Pi 1 Oo) Qf Oo) 

This last condition is what M. Poincare calls the vanishing of a coefficient of 
stability. * It shows that if v 0 and /3 satisfy not only the condition for the Jacobian 
ellipsoid, namely, fjb 1 (v 0 ) (v 0 ) = P, 1 (v 0 ) Q l (v 0 ), but also this equation, we have 

arrived at a figure which belongs at the same time to two series, and there is a 
bifurcation at this point. The form of the figure is found by attributing to e any 
arbitrary but small value. 


SE = — § 


Mp 


Pi 1 to) Qi 1 to) 


§ 6. The Properties of the Successive Coefficients of Stability. 

Corresponding to each harmonic deformation of the ellipsoid, there is a coefficient 
of stability of one of the two forms 

_ IPOo) ©d Oo) __ PdOfJQdOo) 

PJWQJW PJOoWOoV 

These coefficients may be written %■ or K, according to an easily intelligible 
notation. The Jacobian ellipsoid is defined by v 0 , and the question arises as to the 
possibility of the vanishing of the several |£’s as r 0 gradually diminishes from infinity, 
that is to say, as the ellipsoid lengthens. 

An harmonic of the first order merely denotes a shift of the centre of inertia 
along one of the three axes; one of the second order denotes a change of ellipticity 
of the ellipsoid. Since we must keep the centre of inertia at the origin, and since 
the ellipticity is determined by the consideration that the ellipsoid is a Jacobian, 
these harmonics need not be considered, and we may begin with those of the 
third order. 

I shall not attempt to follow M. Poincare in his masterly discussion of the pro¬ 
perties of the coefficients of stability,! but will merely restate in. my own notation 
the principal conclusions at which he has arrived. 

* ‘Acta Math.,’ vol. 7, 1885, p. 321. The factors and 1/2 n + 1 (or 1/21 + 1, if i is the degree of 
the harmonic) which occur in his form of the condition are included in my functions. 

f Sections 10 and 12 of his memoir. I have to thank him for saving me from making a serious mistake 
in this portion of my work, 
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1st. The equation 

Pi 1 ( v ) Qi 1 ( v ) ~ W( v ) <8li s ( v ) or P»* (v) Q i'(v) = 0, (i > 2) 

is not satisfied by any value of v between 1 and infinity, if or P is divisible by 

^ appears from the forms of the functions as given in §4 of 

“ Harmonics ” that the P functions are so divisible. These functions appertain to 
the types EES, OOC, OES, EOO, and therefore the ellipsoid cannot bifurcate into 
deformations of these types. 

2nd. The equation has no solution if IjjJ/ is divisible by (w — l)h We again see 
from § 4 of “Harmonics” that is so divisible if it is of the types OOS, EOS. 
Hence the ellipsoid cannot bifurcate into these types. The only types remaining 
are EEC, OEC. 

3rd. The equation has no solution if any of the roots of (v) = 0 lie outside the 
limits + 1 to — 1. The only 10/ of the types EEC, OEC which has all its roots inside 
the limits + 1 to — 1 is the zonal harmonic for which s — 0. 

Hence the ellipsoid can only bifurcate into a zonal harmonic. 

4th. The equation 

Pi 1 Qi 1 “ = 0 (i > 2) 

must have a solution between 1 and infinity for all values of i. 

It follows from these four propositions that the Jacobian ellipsoid is stable for all 
deformations except the zonal ones, and that as it lengthens it must at successive 
stages bifurcate into each and all the zonal deformations. 

5th. As the ellipsoid lengthens, the first coefficient of stability to vanish is that 
of the third zonal harmonic. This stage is the end of the stability of the Jacobian 
ellipsoids, and there is almost certainly exchange of stability with the pear-shaped 
figure defined by this harmonic. 

6th. It has not been rigorously proved that there is only one solution of the 
equation %i = 0 even in the case where i = 3, but M. Poixcark believes that this 
is almost certainly the case. 

7th. The functions 

U/Ml Si'Wl VW1 WW1 

or > X or > — or > X or > 

PJ«U P/wJ p/MJ p/WJ 

have always the same sign as v increases from v Q to infinity, provided that s and t 
are both greater than zero, and i greater than 2. 

The seventh of the preceding propositions renders it easy to determine the relative 
magnitudes of all the It’s belonging to a single degree i. 

In what follows I may take the symbols |0, <J§. as including also P, Q. 

VOL, CXCVIII. — A. 2 T 
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Now W > = < W as, 

Pd (v 0 ) <©;'* (di) — W (d>) €•/ (Di) < = > 0 ■ 
If we express the ©’s in terms of integrals this becomes 


v i'n r 


The seventh proposition shows that when s and t are greater than Kero, and i is 
greater than 2, all the elements of the integral have the same sign. Hence the 
question is whether. 

W’M ^ __ ^ W("o) 

IfW ^ ^ W«" 


Therefore we have to arrange all the 


(»') 


in 


descending order of magnitude, 


and shall thereby obtain the non-zonal It’s in ascending order, 

I wish first to show that these coefficients may to 'a great extent be sorted by 
considering the inequality. 


I T Oo) ^ ^ l> ‘‘ To) 

IV O') IT(r) 


(s = 1, 2, 3 


i ; t= 1, 2, 3 . . i). 


Suppose, if possible, that whereas, for the ellipsoids defined by /3, v, v () , 

^d(vp) Wbv) , IT (fp) . IT (r 0 ) 
y iVOO TfiW ■ 

Then there must be some value of /3 for which 


ft'WI.‘W = i , W*/W 

for all values of v greater than v 0 . 

It is almost obvious that there is no one value of [?> which renders this equation 
possible; but consider for example the case of s = 2, t = 0. 

Now 

f 3 3 (v) = ~ %>P S ( V ) + P S 2 (*>). % ( V ) = P 3 ( v ) + iSp 2 P 3 2 (v 0 ). 

It we substitute this in the equation we find 

P,>WP,W = P>(,)P S W. 

This can only be satisfied by v — v 0 , and hence the hypothesis is negatived. 
Similarly the assumption of other values of s and t leads to an impossibility. 

Thus we may consider the P functions in place of the ^ functions. 
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Consider the inequality 

IVfa 0 ) 


IV fa) 


fa 


iviLM for 
< iv +1 fa) ’ 1 


1,2. . i - 1. 


If the inequality is determined for any value of v, it is determined for all values. 
Now when v is very large 


p s /\_ 2il 

P® ( V ) 0 i n I «' 


2 1 ill — s! 


JVW-S7 


2 il 


2* Hi 


1 ! 


v. 


Hence our inequality becomes 

(i - s) P,» (* 0 ) > = < Pf 1 (»„). 

This inequality is of the same kind for all values of v 0 . Now P/ (v 0 ) involves the 
factor (v 0 3 — l) |s and P/ + 1 fa 0 ) involves fa 0 3 — l)* (s + 1> . Putting therefore v 0 3 =l + e, 
the left-hand side involves e s and the right + It follows that unless s is equal 
to i the left-hand side is greater than the right; but s is necessarily equal to i — 1 
at greatest. 

Therefore IV fa,) 1 V +1 fa ,) 

P/fa) ^ Pf + 1 fa)' 


Hence K’s with smaller s are less than those with greater s. 

It remains to discriminate between the two sorts of P-functions which occur in 
ellipsoidal harmonic analysis ; that is to say we must determine 

ffiffao) ^ ^ Pf fao) 

Wfa) fa) * 


Since the /3 of “ Harmonics ” is equal to --in the present notation, when fi 

Ju — !C * 

and k are small we have by the formulae of that paper 

w (0 = Pi* (0 + (0 + w + . • •> 


Pi s (0 


fa 3 - lfafa 

“fa 3 - 1) 4 ' 


Pi* {v) + K% +t P^ (v) + K V 2 Pr 2 (0 + . • ■ 


When v is very great and i< very small = P/, so it suffices to determine the 
inequality 

5IT i, v o) V ” V P< ( v o) > 

and this may be considered for any value of v 0 greater than unity. By taking v 0 very 
large and k very small the inequality becomes 

(-= -!)*> = < (v - K )*. 

1 fa = < «. 

2 T 2 


or 
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But k < 1, hence the first sign holds true and 

W_(fo) ^ P / (To) 

whence > IV < K;U 

r riius it follows that for order i 


< K i 1 < mu < K/ . . . . < IV < K/U 

The order of magnitude of these coefficients is therefore completely determined. 

As confirmatory of the correctness of this result it may he mentioned that 1 find 
that when y = 69° 50 7 and k = sin 73° 50', 

Ig 1 = *1765, K 3 1 = *2990, m a 2 = *4467, K,f — *4550, % 6 = *5719, K 3 3 = “5876. 

When y ~~ 75° and k = sin 81° 4 7 (another Jacobian ellipsoid) the numbers run 
•130, *224, *460, *465, *604, *614. 

We see that for the harmonics of higher order the ellipsoid is more stable than it 
was and for those of lower order less stable. 


§ 7. The critical Jacobian Ellipsoid . 

.From a number of preliminary calculations I saw reason to believe that the critical 
ellipsoid would be found within the region comprised between y = 69° 48 7 and 
69° 50', and sin ” 1 k == 73° 52' and 73° 56'. 

If we write 


f(y, sin 1 k) 



(2F _ e) 

1 — k j snr 7 / 


k? sin 7 cos 7(1 + k 2 s in- 7 ) 

1 0 / -1 O * O \X. 

K " (1 — AT Sill 7) 9 


where the amplitudes of E and E are y and their moduli k, the existence of the 
Jacobian ellipsoid is determined by 


f(y, sin 1 k) — 0 .* 


The coefficient of stability is 


% (y, sin 1 K ) — 1 — 


IPs Oo) @3 M 
P/KlQJO'o)' 


The formulae for computing 1& 3 are given in § 4. 

The values of E and F are from Legendre’s tables. 


* See ( Roy. Soc. Proc., 5 voL 41, p. 323 x where the formula is reduced to a form convenient for 
computation. 
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Now I find 

/(69° 48', 73° 52') = + -000191 ; /(69° 50', 73° 52') = + -001319. 

/(69° 48', 73° 56') = - -001186; /(69° 50', 73° 56') = - -000031. 

% (69° 48', 73° 52') = + ’001058 ; E 3 (69° 50', 73° 52') = - '000885. 

3stg (69° 48', 73° 56') = + ‘000655 ; lt 3 (69° 50', 73° 56') = - ’000765. 

By interpolation we get the following results :— 

The Jacobian ellipsoid is given by 

(y — 69° 48') — -59642 (sin" 1 k - 73° 52') + *33091 = 0. 

The vanishing of the coefficient of stability is given by 

(y - 69° 48') + -041625 (sin" 1 k - 73° 52') - 1*0890 — 0 . 

In these equations the minute of arc is the unit. 

Solving them I find 

y = 69° 48'-997 = 62° 49'*0, 
sin -1 * = 73° 54'-225 = 73° 54'*2. 

With these values I find that the three axes a, b, c, where abc = a y are 

“ — -650659, 

•f- — -814975, 

-- = 1-885827. 

a 

The last place of decimals in these is certainly doubtful. 

The formula for o>' 1 is given in (35). 

Now IF = | Po ( v o) ©o ( v o)> li — CK siu r> Po ( v o) ©o ( v o) = K ^- 

Then since a = c cos y, b = cA, 

M 2 __ 2F A cot 7 - 
2rrp 1 + A 3 - 1T5 =rn5^ * 

In this formula, F, y, A must correspond with values interpolated amongst those 
used in obtaining the solution. 

From this I find 

= -1419990 = -14200. 

27 Tp 
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In the paper on the Jacobian ellipsoid referred to above the moment of momentum 
is tabulated by means of p, where the moment of momentum is (firp)ia 5 p. The 
formula for p is given in (25) of that paper, and, modified to suit the present 
notation, is 

p = |(4eo Sr )-i(I+A s )(£r (( )'- 

For the critical ellipsoid 1 find /x = ‘889570. 

The following table gives the numerical values for a number of Jacobian ellipsoids, 
beginning with the initial one and terminating just beyond instability. The last line 
gives the corresponding values for the critical ellipsoid. 


Jacobi’s Ellipsoids.* 


7* 

sin~ 1 K, 

cos “ a A. 

a/a. 

If/ a. 

c/a. 

« 3 /27rp. 

/u 

54 21 27 . . 

0 0 

0 0 

•6977 

IT 972 

1T972 

•18712 

30375 

55 .... . 

17f 

14-1 

•697 

1*179 

1*216 

T8706 

*304 i 

57 .... . 

34f 

28| 

*696 

1*123 

1*279 

*186 

: *306 

60 .... . 

49 7 

40 54 

*6916 

1*0454 

1*3831 

*1812 

*3134 

65 , 

64 19 

54 46 

•6765 

*9235 

1*6007 

*1659 

*3407 

70 .... . 

74 12 

64 43 

•6494 

•8111 

j 1*899 

j *1409 

*3920 

69 49 

73 54 

64 24 

65066 

*81498 

1*88583 

*14200 

*38957 


* I have been criticised with respect to my paper on Jacobi's ellipsoid, from which these results are 
extracted, by M. S. Kruger (Nieuw Archief voor Wiskunde, Tweecle Keeks, Derde Deel and 
‘ Ellipsoidale Evenwichtsvormen/ &c., Thesis for Degree of Doctor, Leiden, J, W. van Leeuwen, 
Hoogewoerd 89, 1896), because I wrote it in ignorance of certain previous work, especially of a paper 
by Plana (‘ Ast. Nadir.,’ 36, n. 851, c. 169). But I cannot but congratulate myself on my ignorance, 
since it appears that Plana gave a number of numerical results which were wholly wrong, A knowledge 
of that paper would no doubt have caused me much further trouble. 

My paper gives a number of solutions of the problem which I believe to be correct. Unfortunately 
the methods of the paper are clumsy, and there are several mistakes. The formula for to 2 used in this 
present paper, is much better than that given there. 

The complicated formula on p. 325 is susceptible of reduction to a simple form, for on substituting for 
y its approximate form (i) we have simply 

y - 8 = \ k 2 sin S cos 6, 

where 6 - 54° 2V 27". 

The final numerical result was, however, nearly right, for I now find 

sin = io*9266S2i gin (y - 6), 

whereas I had *9266528. The sin sc is the same as the k used here. 

The formula at the top of p. 326 which is reproduced as (22) on p, 828 is, I think, illusory, for if in the 
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In order to determine the question as to whether or not it is possible that Ms — 0 
should have another solution than that found in the next section, I have computed 
the value of this coefficient for the Jacobian ellipsoid y — 75°, k = sin 81° 4 '' 4 , and 
find it to be —6‘627. From the manner in which the numbers in the computation 
present themselves, it is obvious that for more elongated ellipsoids Ms will always 
remain negative, and will become numerically greater. I have therefore not thought 
it necessary to seek for an algebraic proof that there is no second root of the 
equation. 

Very long Jacobian ellipsoids tend to become figures of revolution, and the coeffi¬ 
cients of stability tend to assume the forms 

, _ PiO) Q»(r) 

Px'O'WW* 

The forms of these functions are well known, and I think that fair approximations 
to the incidences of the successive figures of bifurcation might be derived from the 
vanishing of this expression. 

For example 

Pi 1 M Q. 1 M = *[■'-(** — i) log (— i)'' 

FtMQ.M = wi (sue - + 3 > log fiAt)'- (21v ! - 11)(»S»*— 30S + 3) . 

I have not, however, attempted to solve the equation found by equating these two 
expressions to one another. 

Even when i = 3 and y = 69° 49' (the critical Jacobian) this rough approximation 
makes the coefficient of stability very small, but it is to be admitted that P^ Q x x and 
P 3 Q 3 differ very sensibly from P x l (v) Q x 3 ( v ) and (z^) (©3 ( v ), although in such a 
way that the errors compensate one another. 

first term we put y = 8 + \ k 1 sin 8 cos 8 (as is clearly allowable in approximation) the term with eoeffi- 
cient k~ or sin 2 a disappears. This shows that it was necessary to proceed in the approximation as far as 
k 4 in order to obtain a result. 

The methods of approximation adopted on pp. 326-7 are of doubtful propriety, but will, I think, lead to 
nearly correct results. There is, however, a mistake towards the bottom of p. 327 which runs on to the 
end. M. Kruger correctly points out that the second line of formula (24) p. 329 should run 

-i cos 2 a ~J~ loge cot (|-7r - |y), (- 1 / 4 - 3 tan 2 y + tan 4 y) ~ - 1 / - - 2 / tan 2 y -1- tan 4 y . 

u sm y ___ 

Lastly, on p. 335, line 13, for 0 = 0*3573, read 0 = 0*5379; and on p. 336, line 7, for 1*3573, read 

1*5379; and for - = 1*696, read- = 4*65. 

a a 
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§ 8. The pear-shaped Figure of Equilibrium. 

By (21) the normal displacement Sn for the third zonal harmonic deformation may 
be written 

o s[gV-g»* »-(3 - 

«Y* (1 - ? 3 Rin 3 7) (a?/coB 4 y + f /A\ + **)* ’ 

subject to the condition 


--'■V + ^ + = C* 

cos" 7 A" 


The expression has been arranged so that when x = y — 0, z = o, we have 8n = e. 
Hence + e and — e are the normal displacements at the stalk and blunt end of the 
pear respectively. 

In the section y — 0, this may be written 

_ 6 cos v — c ¥0 
dn ~ <f- • - * siii® yf ‘ 


The nodal points are given by = rt 9 = ± ’758056. 


4 - 5 cf 


In the section x = 0, since k = p 3 -—— , it may be written 


8??, 


4 f 

A (4 — 5(f) z(«% 3 ~cY) 

f- .* cV(e 3 - «W ! sin- 7 ; J 


• (i 

The nodal points are given by j = fit - = ± ’788986. 

The section z — 0 is obviously another nodal line for all sections. 

By means of these formulae it is easy to compute the normal displacements from the 
surface of the critical Jacobian. 

The figure opposite showing the three sections * = 0, y = 0, z = 0, is drawn 
from these formulae, the dotted line being the critical Jacobian and the firm line the 
pear. The scale of the normal displacements is, of course, arbitrary. 

Comparison with M. Poincare’s sketch shows that the figure is considerably longer 
than he supposed. 

In this first approximation the positions of the nodal lines are independent of the 
magnitude of e, and they lie so near the ends that it is impossible to construct an 
exaggerated figure, for if we do so the blunt end acquires a dimple, which is absurd. 
It might have been hoped that such an exaggeration would afford us some idea of the 
mode of development of the pear. 
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M. Schwarzschild has remarked* that it is not absolutely certain that the 
principle of exchange of stability holds with reference to this figure, and that we 
cannot feel absolutely certain that the pear is stable unless we can prove that the 
moment of momentum is greater than in the critical Jacobian. 

With reference to this objection, M. Poincare writes to me as follows :— 

“ Faisons croltre le moment de rotation, que j’appellerai M. Deux hypotheses sont 
possibles. 

“ Ou hien pour M < M 0 (the moment of momentum of the Jacobian), nous aurons 
une seule figure, stable, a savoir l’ellipsoide de Jacobi, et pour M > M 0 trois figures, 
une instable, l’ellipsoide, et deux stables (d’ailleurs egales entre elles), les deux figures 
pyriformes. 

“ Ou bien pour M < M 0 , nous aurons trois figures d’equilibre, deux pyriformes 
instables, une stable, l’ellipsoide, et pour M < M 0 une seule figure instable, l’ellipsoide— 
auquel cas la masse fluide devrait se dissoudre par un cataclysme subit. 

“Ilya done a verifier si pour les figures pyriformes, M > ou < M 0 .” 

It seems very improbable that the latter can be the case; but this opinion is not a 
proof. 

Since co 2 is stationary for the initial pear, a small change in the angular velocity 
will certainly produce a great change in the figure of the pear. If this investigation 
has, in fact, its counterpart in the genesis of satellites and planets, it seems clear that 
the birth of a new body, although not cataclysmal, is rapid. 


§ 9. Summary. 

It is possible by the methods explained in my previous paper on “ Harmonics” to 
form rigorous expressions for the ellipsoidal harmonics of the third degree. Accordingly 
in § 1 I proceed to form those functions. In § 2 the notation is changed with a view 
to convenience in subsequent work, and for the sake of completeness the harmonics 
of the first and second degrees are also given. In § 3 the corresponding solid 
harmonics are expressed in rectangular co-ordinates x, y, z. In § 4 I find the 
Q-functions, the harmonic functions of the second kind, and express the results in 
terms of the elliptic integrals E and F. It appears that both the P- and Q-functions 
of the third degree of harmonics occur in three pairs which have the same algebraic 
forms, and that in each pair one of them only differs from the other in the value of a 
certain parameter. There is, lastly, a seventh function which stands by itself; this 
last corresponds to the solid harmonic xyz. 

In § 5 the equations for Jacobi’s ellipsoid are determined by the consideration that 
the energy must be stationary, and the superficial value of gravity is found in terms 
of the appropriate P- and Q-functions. I then proceed to find the additional terms 

* “ Die Poincaresche Theorie des Gleichgewichts,” ‘ Aimalen der K. Sternwarte, Miinchen,’ Bd. III. 
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in the energy when the mass of fluid is subject to an ellipsoidal harmonic deformation. 
This section is a paraphrase of M. Poincare’s work, but the notation and manner of 
presentation are somewhat different. The additional terms in the energy are shown 
to involve a certain coefficient, which is called by M. Poincare a coefficient of 
stability. It is clear that when any coefficient vanishes we are at a point of bifurca¬ 
tion, and the particular Jacobian ellipsoid for which it vanishes is also a member of 
another series of figures of equilibrium. 

In § 6 the principal properties of these coefficients, as established by M. Poincare, 
are enumerated. He has shown that the ellipsoid can bifurcate only into figures 
defined by zonal harmonics ; that it must do so for all degrees, and that the first 
bifurcation occurs with the third zonal harmonic. The order of magnitude of the 
coefficients of the several orders and of the same degree is determined. A numerical 
result seems to indicate that as the ellipsoid lengthens, it becomes more stable as 
regards deformations of the third degree and of higher orders, and less stable as 
regards the lower orders of the same degree. 

In § 7 the numerical solution of the vanishing of the coefficient corresponding to 
the third zonal harmonic is found, and it is shown that the critical ellipsoid has its 
three axes proportional to '65066, '81498, 1'88583, and that the square of the angular 

velocity is given by = '11200. A short table is also given showing the march 

of the axes of the Jacobian ellipsoids from their beginning on to instability at this 
critical stage. The nature of the formula for the third zonal coefficient of stability 
seems to show that it can only vanish once—a point which it appears that 
M. Poincare found himself unable to prove rigorously. 

A suggestion is made for the approximate determination of the bifurcations into 
the successive zonal deformations, but no numerical resxdts are given. 

In § 8 the nature of the pear-shaped figure is determined numerically, arid the 
reader may refer to the figure above, where it is delineated. It will be seen to be 
longer than was shown in M. Poincare’s conjectural sketch. 

If, as M. Poincare suggests, the bifurcation into the pear-shaped body leads 
onward stably and continuously to a planet attended by a satellite, the bifurcation into 
the fourth zonal harmonic probably leads unstably to a planet with a satellite on 
each side, that into the fifth to a planet with two satellites on one side and one on the 
other, and so on. 

The pear-shaped bodies are almost certainly stable, but a rigorous and conclusive 
proof is wanting until the angular velocity and moment of momentum corresponding 
to a given pear are determined. To do this further approximation is needed. 


2 u 2 



